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In the final paper of tfiis series, we extend our results on magnification invariants to the infinite 
family of An {n > 2),D„ (n > 4), Eg, E7, Eg caustic singularities. We prove that for families of 
general mappings between planes exhibiting any caustic singularity of the An {n > 2) , Z)„ (n > 
i), Ee, E7, Eg family, and for a point in the target space lying anywhere in the region giving rise 
to the maximum number of lensed images (real pre-images), the total signed magnification of the 
lensed images will always sum to zero. The proof is algebraic in nature and relies on the Euler trace 
formula. 
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I. INTRODUCTION 

^ I In papers I and II of this series (Aazami & Fetters 2009 we established a universal maenification theorem 

for all higher-order caustics up to codimension 5. It was shown that to each such caustic singularity is associated a 
magnification sum rule of the form 
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^971, =0. 



PsJ All of these bifurcation sets or big caustics occur in a n-parameter space. Slices of the big caustics give rise to caustic 
K>" metamorphoses that occur in gravitational lensing (e.g., Blandford 1990 J^, Schneider, Ehlers, & Falco 1992 jT3| , 
^ Fetters 1993 [11], and Fetters et al. 2001 Chaps. 7, 9]). It was discussed in Faper I, using the hyperbolic umbilic 
in particular, how the above magnification relations may be used for substructure studies of four-image lens galaxies. 

We now extend these results to all higher-order caustic singularities beyond codimension 5. These are classified 
according to Arnold's A,D,E classification of Lagrangian map-germs (Arnold 1973 [3]). Thus, we show that for 
families of general mappings between planes exhibiting any caustic singularity with corresponding Coxeter-Dynkin 
diagram of type An (n > 2), D„ (n > 4), Eq, E^, Eg, and for a point anywhere in the region of the target space giving 
I rise to the maximum number of lensed images (real pre-images), the total signed magnification is identically zero. Our 
proof is algebraic and relies on the Euler trace formula. Finally, we emphasize that the magnification sum relations 
are geometric invariants, being the reciprocals of Gaussian curvatures at critical points. 

The outline of this paper is as follows. In Section |lT] we give a brief overview of the A, D, E family of caustic 
singularities. In Section [IIII we state our main theorem. The proof itself is presented in Appendix [Xj 



II. HIGHER-ORDER CAUSTICS OF THE A, D, E FAMILY 



We inherit the notation and terminology of Faper II. To that end, consider a smooth n-parameter family i^c,s(x) of 
functions on an open subset of that induces a smooth (n — 2)-parameter family of mappings fc(x) between planes 
{n > 2). One uses F^^s to construct a Lagrangian subnianifold that is projected into the space {c,s} = x R^. 

The caustics of fc will then be the critical values of the projection (e.g., Golubitsky & Guillemin 1973 [9], Majthay 
1985 [l3l, Castrigiano & Hayes 1993 0, and [l2l . pp. 276-86]). These projections are called Lagrangian maps, and 
they are differentiably equivalent to fc. 

Arnold classified all stable simple Lagrangian map-germs of n-dimensional Lagrangian submanifolds by their gen- 
erating family i^c,s ([3], Arnold, Gusein-Zade, & Varchenko I 1985 ^ p. 330-31], and [l^, p. 282]). In the process 
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A„ (n > 2) 


Fc,six,y) = ±x"+^ ±1/2 +c„_ix"-^ H hC3X^ + S2X^ - s^x ± S2y 

germ unfolding terms 

fc(x,y) = (±(n+l)x" + (n-l)c„_ix"-2 + --- + 3c3x2=F4|/x , =F 2?/) 


D„ {n > 4) 


Fc,s(x, y) = x^j/ ± y""^ + c„_2j/"~2 H h £2?/^ - S2y -six 

germ unfolding terms 

fc(x, y) = (2xy , ± {n - + (n - 2)c„_2y"-^ + ■ • • + 2c2y) 


Ee 


Fc,b{x, y) = x^ ±y'^ + C3ZJ/2 + £22/^ + cixy — S2y — six 

germ unfolding terms 

fc{x, y) = {3x^ + c-iy'^ + ciy , ± Ay'^ + 2c3xy + 2c2y + cix) 


Er 


Fc,s{x,y) = x^ + xy^ + C4j/* + c^y^ + C2y^ + cixy - S2y — six 

germ unfolding terms 

fc{x, y) = (3a;2 + + ciy , Sxj;^ + Ac^y'^ + Scgy^ + 2c2j/ + cix) 


Es 


Fc,s{x, y) = + y^ +C5xy^ + C4xy'^ + csy^ + C2J/2 + cix-y - S2y - six 

germ unfolding terms 

fc{x, {/) = (3x2 _|_ ^^j^3 (,^j,2 ^ ^ 5j^4 ^ 3cgxj/2 + 2c4Xj/ + Scaj/^ + 2c2j/ + cix) 



TABLE I: For each type of Coxeter-Dynkin diagram listed, indexed by n, the second column shows the corresponding universal 
local forms of the smooth (n — l)-parameter family of general functions Fc,s, along with their (n — 3)-parameter family of 
induced general maps fc between planes. This classification is due to Arnold 1973 



he found a connection between his classification and the Coxeter-Dynkin diagrams of the simple Lie algebras of types 
An {n > 2),Dn [n > 4), E^, Et, Eg. This classification is shown in Table HI The singularities in Paper II arose as 
follows: A2 (fold), A3 (cusp), (elliptic umbilic), (hyperbolic umbilic), A4 (swallowtail), Ar-, (butterfiy), D5 
(parabolic umbilic), Aq (wigwam), Eq (symbolic umbilic), (2"^^ elliptic umbilic), and (2""^ hyperbolic umbilic). 

For the shown in TablelH call x G a lensed image (or a real pre-image) of the target point s G if fc(x) = s 
(in particular, our lensed images are always in M^, not C^). A point x^ G is a lensed image of the target point 
s G if and only if x^ is a critical point of i^c,s (relative to a gradient in x). Next, we define the magnification 
3Jl(xj; s) at a critical point Xi of the family i^c,s by the reciprocal of the Gaussian curvature at the point (xi, i^c.s(xi)) 
in the graph of Fc^s- 

»t(x, ; s) = , ^„ , • (1) 

Gauss(xi,Fc,s(xi)) 

This makes it clear that the magnification invariants established in our theorem are geometric invariants. In addition, 
since Gauss(xi, i^c,s(xi)) — det(Hess Fc,s)(xi), and since each fc in Table |T] satisfies det(Jacfc) — det(Hess Fc,s), we 
can also express the magnification in terms of fc: 

mt(x,;s) ^ 



det(Jac fc)(xi) 



where x^ is a lensed image of s under fc. We now proceed to our main theorem. 
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III. STATEMENT OF MAIN THEOREM 



Theorem 1. For any of the universal, smooth [n — \) -parameter family of general functions i^c.s (or induced general 
mappings f,,) in Table\^ and for any non-caustic point s (light source position) in the indicated region, the following 
results hold for the magnification 931^ = ^^(xi ; s) ; 

1. An {n > 2) obeys the magnification relation in the n-image region: 'Y^^^i = 0; 

2. Dn (n > 4) obeys the magnification relation in the n-image region: X]r=i ^« ~ ^' 

3. obeys the magnification relation in the six-image region: X]i=i ~ 

4- Ej obeys the magnification relation in the seven-image region: X]I=i — 0, 

5. Es obeys the magnification relation in the eight-image region: Y^^=i ^« ~ ^■ 

The magnification relations for A2 (fold) and As (cusp) are known [gI. [l^ - fTsj . Those for (elliptic umbilic), 
(hyperbolic umbilic), A4 (swallowtail), (butterfly), (parabolic umbilic), Aq (wigwain), Eq (symbolic umbilic), 
Dq (2"*^ elliptic umbilic) , and (2"^^ hyperbolic umbilic) were discovered recently in [l|, |2[ . 

Remarks. First, the results of Theorem [1] actually apply even when the non-caustic point s is not in the maximum 
number of real pre-images region. However, pre-images from will appear, which are unphysical in gravitational 
lensing. Second, it is important to point out that Theorem [1] is not a direct consequence of the Euler- Jacobi formula, 
of multi-dimensional residue integral methods, or of Lefschetz fixed point theory because some of the singularities 
have fixed points at infinity that must be treated invidually. Third, we point out that for n > 6 there are Lagrangian 
maps that cannot be approximated by stable Lagrangian map-germs jSi] . 



The paper presented a theorem about the magnification of lensed images for all caustic singularities appearing in 
the infinite family of An (n > 2),£)„ (n > A) , Eq , E-j , Eg, caustic singularities. We proved that for families of general 
mappings between planes locally exhibiting any caustic singularity of the An (n > 2),Z?„ (n > 4), Eq, Ey, Eg family, 
and for a target point lying anywhere in the region giving rise to the maximum number of lensed images (real pre- 
images), the total signed magnification of the lensed images will always sum to zero. The signed magnifications are 
geometric invariants as they are Gaussian curvatures at critical points. The proof was algebraic in nature and made 
use of the Euler trace formula. Our result goes beyond previous work that considered singularities up to codimension 
five. 



AOP acknowledges the support of NSF Grant DMS-0707003. Part of this work was conducted at the Fetters Research 
Institute, Belize. 



We summarize key elements of our algebraic method; see Faper II for a detailed presentation. Consider any polynomial 
(p{x) = Onx" + ■ ■ ■ + aix-\-ao S C[x] with distinct roots xi, . . . ,Xn and any rational function h{x) 6 i?, where R C C{x) 
is the subring of rational functions defined at the roots of (p{x). In Faper II we showed algebraically that for any 
rational function h{x) £ R, the following holds: 



IV. CONCLUSION 
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Appendix A: Proof of the Main Theorem 



1. 



Overview of the Method and the Euler Trace Formula 




(Euler Trace Formula) 



(Al) 
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where 6„_i is the (n — l)st coefficient of the unique polynomial representative r{x) of degree less than n in the coset 
(p'{x) h(x) G R/{tf{x)). (An alternate proof of the Euler trace formula using residues can be found in Dalai & Rabin 
2001 Q.) We employ the Euler trace formula as follows. For any caustic singularity in Table HI we shall realize the 
lensed images of its corresponding family of mappings fc as solutions of a polynomial in one variable, obtained by 
eliminating one of the pre-image coordinates, say y. Denote this polynomial by 

1^9(2;) = a„a;" + • ■ • + aia; + ao G C[a;]. 

Generically, we can assume that the roots of Lp{x) are distinct, an assumption made throughout the paper. We would 
then be able to express the magnification S[Jl(a;, y; s) at a general pre-image point {x, y) as a function of one variable, 
in this case x, so that 

mix, yixy,s) = det(Jacfe)(a;,y(a;)) ^ det(Jac fc)(x) ^ ^^""^ ' 

where the explicit notational dependence on s is dropped for simplicity. Since we shall consider only non-caustic target 
points s giving rise to lensed images {xi,yi{xi)), we have det(Jac fc)(xi) 7^ 0. We thus know that for the singularities 
in TablelH the rational function 9Jl(x) is defined at the roots of f(x), i.e., dJl{x) € R. Now, denote by Tn(a;) the unique 
polynomial representative of degree less than n in the coset Lp'{x)dJl{x) G R/{ip{x)), and let 6„_i be its (n — l)st 
coefficient. In the notation used above, we have h{x) = 2n(x) and r{x) = m{x). Euler's trace formula (Corollarv lAip 
then tells us immediately that the total signed magnification satisfies 

(A2) 



It therefore remains to determine the coefficient bn-i for each caustic singularity in Table HI 



2. Magnification Sum Rule for Type An 

We begin with type An, n > 2. Since the cases 2 < n < 6 were treated in Paper II, we will consider n > 7 here. The 
(n — l)-parameter family of general functions F^" is given in Q by 

F^" {x, y) = ±a;"+^ ± + c„_ia:""i + c„_2a;""^ H h cgx^ 4- C2X^ + cix . (A3) 

V ' V ' 

germ unfolding terms 

To convert this into the form shown in Table HI we use the following coordinate transformation on the domain 
{(x,y)}=R2: 

{x,y)^ (^x,y+^) . (A4) 

This transforms eqn. (jA3p to 

F^si^^y) = ±a;"+^ ±y^ + Cn-ix"-'^ + c„_2a;""^ + ■ ■ • + cs.x'"^ + S2x'^ - Six ± S2y , (A5) 

where ci = — si and C2 = 82- The parameters si, S2 are to be interpreted in the context of gravitational lensing as the 
rectangular coordinates on the source plane 5 = . Note that we omitted the constant term from eqn. (|A5[) since it 
will not affect any of our results below. Note also that 

det (Hess F^" ) = det (Hess F^^ ) , 

so that the magnification (as defined in eqn. ([1])) is unaltered. We will work with the form of F^^ in eqn. (jA5l) . The 
corresponding {n — 3)-parameter family of general mappings f^" : — > R^ is 

i^"{x,y) = {±{n + l)x'' + {n - l^n-ix"'^ + {n - 2)cn-2x''-^ + ■ ■ ■ + 3c3X^ T 'iyx , T^y) = (si,S2). 

Here s = (si, S2) is a non-caustic target point lying in the region with the maximum number of lensed images. Since 
S2 = ±2y, we can eliminate y to obtain a polynomial in the variable x: 

ipA„ (x) = ±{n + l)a;" + {n - l)c„_ix"~2 + (n ~ 2)cn-2x"^^ H h Scsx'^ + 2s2X - si , (A6) 
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whose n roots are the x-coordinates of the lensed images of s. The Jacobian determinant of f^" expressed in the 
single variable x is 



det (Jacf^^") {x) = t2 [±n{n + l)x"-^ + {n ~ 2){n - l)c„_ix""^ + (n - 3)(n - 2)cn-2x"^^ + ■■■+ Gc^x + 2s2] 



(A7) 



A comparison of eqns. (|A6P and (|A7P then shows that 



We thus have 



±2(p'a^{x) = det (Jacf^") (a;) 



^'^^{x)m{x) = ±-- 



m{x) 



Thus the unique polynomial representative of the coset (x)93t(a;) is the polynomial m{x) = ±1/2, whose {n — l)st 
coefficient is 6„_i = for all n >7. Euler's trace formula in the form of eqn. (jA2[) then tells us that the total signed 
magnification is 

n 

J2^^=0, {A„,n>2). 



3. Magnification Sum Rules for Type D„ 

For type > 4, the cases 4 < n < 6 were treated in Paper II, so we will consider n > 7 here. The corresponding 

(n — 3)-parameter family of induced general maps fc " : — > 

shown in Table ID 

" (x, y) = (2xy , x^ ± {n ~ l)y"-2 ^ „ 2)c„_22/"-=' + • ■ • + (n - z)c„_,2/"-(^+i) + • • • + 2c2y) = (si, S2). (A8) 

Once again the point s = (si, S2) is a non-caustic target point lying in the region with the maximum number of lensed 
images. This time, however, we eliminate x to obtain a polynomial in the variable y: 

ip^± (y) = ±4(n - 1)2/" + 4(n - 2)c„_2j/""' + • • • + 4(n - i)c„_,y"-(*-i) + ■ • • + 8c2y3 - 4s2y^ + , 

whose n roots are the y-coordinates of the n lensed images Xi of s. The derivative of Lp^±{y) is 

ip'i (y) = ±4n(n - l)y"-i + 4(n - l)(n - 2)c„_22/"-' + ■ • ■ + 4(n - - l))(n - z)c„_,y"-^ + ■ • ■ + 24c22/' " 8^22/ , (A9) 



while the Jacobian determinant of fc " is 
det ^Jac ftP" ^ (x, y) = det 



2y 2a; 
2a; ±(n - 2)(n - l)y"-3 + (n - 3)(n - 2)c„_22/""^ + 
±2(n - 2)(n - l)y"-2 ^ 2(n - 3)(n - 2)c„_2y""^ + • • ■ 
• • • + 2(n - (^ + l))(7i - ^)c„_,y"-(*+l) + • ■ • + 4c2y - 4a;^ 



2C2 



We can use eqn. (|A8[) to eliminate a; as follows: 

= ±2(n - 2)(n - l)y"-2 ^ 2(n - 3)(n - 2)c„_2y""'' + ■ • ■ + 2(n - + \))(n ~ z)c„_,y"-(*+i) + ■ • • + 4c2y 
+4 (±(n - l)y"-2 ^ _ 2)c„_22/"-' + ••• + («- i)c„„,y"-(^+i) + • ■ • + 2c2y - 



: —x^ (by oqn. ( lASH 



= ±2n{n ~ l)y"-^ + 2(n - \){n - 2)c„_2y"-^ + • • ■ + 2(n - (i - l))(n - i)c„_,y"~('+i) + • ■ • + 12c2y - 4s2 



det (jacff") (y) = an(y)- 
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A comparison with eqn. (jA9[) then shows that 

The unique polynomial representative of the coset '/'^i iy)^{y) is therefore the polynomial m{y) = 2y, whose (ri — l)st 
coefficient is = for all n > 7. Eqn. (IA2p then tells us that the total signed magnification is 

n 

^971, = , p„, n>4). 



4. Magnification Sum Rules for Types En 

The case Eq corresponds to the symbolic umbilic, whose magnification sum rule was proved in Paper 11. For type iJy, 
Table |T] gives the corresponding 4-parameter family of induced general maps f^'' : — > : 

ffHa;,y) = {3x^ + + ciy , 3xy^ + Aay^ + Scsy^ + 2c2y + cix) = (si,S2). (AlO) 

Once again, the point s ~ (si,S2) is a non-caustic point lying in the region with the maximum number of lensed 
images. We eliminate x to obtain a polynomial in the variable y: 

ifEriy) = + 48cl/ + (15ci + 72C3C4) y^ + {274 + 48C2C4 - 9si) y^ + (7c? + 36C2C3 - 24C4S2) y^ 
+ {I2cl - 6cisi - I8C3S2) 2/^ + {c\ ~ I2C2S2) y - clsi + 3s^ , 

whose 7 roots are the y-coordinates of the 7 lensed images of s. The derivative of (pEriy) is 
ip'j^^iy) = 63/ + 288c^y5 -|_ (jSc^ + 36OC3C4) y^ + (l08c^ + 192c2C4 - 36si) y^ 

+ (21c? + IO8C2C3 - 72C4S2) y^ + (24c^ - 12cisi - 36C3S2) y + 4 - 12c2S2 , (All) 

while the Jacobian determinant of ff'" is 



det (Jacff^) {x,y) = det 



6x 3y^ + ci 

Sy^ + ci 6xy + 12c4y^ + 603?/ + 2c2 



— 36x^y + 72cixy^ + SGc^xy + 12c2X — 9y'^ — 6ciy^ — c{. 
To convert this into a function in the single variable y, we use eqn. (|A10[) twice, as follows: 

= i&x^y + 72cixy^ + 36c3a;y + 12c2a; — Ot/*^ — Qiy^ — c\ 

si-ciy~y^\ , ^^^_.2 , o«„ „, , io„ \ /^*2 - 4c4y3^-3c3y2 - 2c2j/^ „4 2 „i 

3j/2 + ci 



362/ ( ) + (72c4y^ + i^c^y + 12c2) ( , , ) -9y" - 6ciz/ - c^ 



= a:2 (by eqn. llAlOl 'l =x (by eqn. lAlOt ) 

= (ci + iy^y^ [ - c? + 12c2S2 + ("24c^ + 12ciSi + 36C3S2) y + (-21c? - IO8C2C3 + 72C4S2) y^ 

+ (-108c^ - 192C2C4 + 36si) y^ - (75ci + 36OC3C4) y'' - 2%%c\y^ - 63y^] 
= det(Jacff^)(?;) = 9Jt(2/)-^ (A12) 

A comparison of eqns. (jAlip and (|A12|) then shows that 

(p's^(y)an(y) = -(ci + 3y2) , 



and so the unique polynomial representative of the coset (/3^„(?/)9Jt(j/) is the polynomial m(y) = — (ci + 3?;^), whose 
6th coefficient is 65 = 0. We conclude via eqn. (|A2|) that the total signed magnification is 

7 

i=l 
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For type Es , Table U gives the corresponding S-parameter family of induced general maps f^f : K.^ — !• : 

i^'{x,y) = {3x^ + C51/ + C41/ + ciy , + Sc^xy^ + 2c4xy + Scsy^ + 2c2y + cix) = (si,S2). (A13) 

As usual, the point s = (si, S2) is a non-caustic target point lying in the region with the maximum number of lensed 
images. We eliminate x to obtain a polynomial in the variable y, 

(fiEsiy) = 75?/® + 9c5?/^ + (90c3 + 21c4cl) + (60c2 + IGc^cg + Ibcicj) y^ + (27c^ + 4c| + 22C1C4C5 - 9c^sl - SOsz) 
+ (36C2C3 + 8C1C4 + 7clc5 - 12C4C5S1) j/^ + (l2c2 + 5clc4 ~ 4C4S1 — 6C1C5S1 - I8C3S2) 
+ (cj - 4ciC4Si - I2C2S2) y - cjsi + 3sl , 

whose 8 roots are the y-coordinates of the 8 lensed images of s. The derivative of fEsiu) is 
VEsiy) = 600/ + eScsy*^ + (540c3 + 126c4c2) y5 + (300c2 + 80c^C5 + 75cic2) y4 

+ (I08c^ + 16cl + 88C1C4C5 - 36c^si - 12OS2) + (IO8C2C3 + 24cic^ + 21c?c5 - 36C4C5S1) 



24c2 + 10c]^C4 — 8C4S1 — I2C1C5S1 — 36C3S2) y + Ci ~ 4C1C4S1 — I2C2S2 



(A14) 



and the Jacobian determinant of f^** is 

det (Jacff**) {x,y) = det 



6x Sc^y^ + 2c4y + ci 

3c5i/^ + 2c4y + ci 20y"^ + 6c5xy + 2cax + 6031/ + 2c2 

+ 12c2a; + 12c4X^ — 4ciC4?/ + 36c3a;y + 36cc,x^y — ^c^y^ 
6ciC5y^ - 12c4C5y^ + 120a;y^ - 9c^y*. 



Similar to the case -E7 above, we convert this to a function in the single variable y with the aid of eqn. (jAT3j: 
= —cf + \2c2X + 12c4X^ — 4ciC4?/ + i(!>czxy + 3&c^x^y — 'ic^y^ 



-6ciC5y^ - 12c4C5y^ + 120a;y'^ - 9c^y^ 



12c2 + 36c3y + 120y^ 



S2 - 5y^ ~ 3c3y'^ - 2c2y 
3c5y2 + 2c4y + ci 



(12c4 + 36c5y) 



si - ciy ~ Ciy'^ - C5y^ 



= X (by oqn. IIA13t 'l 

-c\ - 4ciC4y - 4c4y^ - 6ciC5?/^ - 12c4C5y^ - Qcly'^ 



-x^ (by oqn. iIATsI 'I 



= (ci + 2c4y + 3c5y^)"^ [ - + 4ciC4Si + 12c2.S2 + (-240^ - 10c^c4 + 8c|si + YIc^c^sy + 36C3S2) y 
+ (-IO8C2C3 - 24cic^ - 21c?c5 + 36C4C5S1) y^ + (-108c^ - 16c| - 88C1C4C5 + 36c^si + I2OS2) y^ 
+ (-300C2 - 80c|c5 - 75cic^) y'' - (540c3 + 126040^) y^ - 63c^y^ - 600y^] 

^ det (Jacff«) (y) = 9Jl(y)"^ (A15) 

A comparison of eqns. (|A14I) and (jA15[) thus shows that 

V'k(y)97J(2/) - -(ci + 2c4y + 3c5y') , 



and so the unique polynomial representative of the coset Lp'-^^ {y)'iffl{y) is the polynomial rn(y) = — (ci + 2c4y + 3c5y^), 
whose 7th coefficient is 67 = 0. Euler's trace formula tells us once again that 

8 

^OH, = 0, (Es). 



This completes the proof of the magnification relations for all higher-order caustic singularities appearing in the 
infinite A„ {n > 2), £)„ {n>4),EG,E7, Es family. □ 
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